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ABSTRACT 

This  study/  assuming  stationary  Gaussian  turbulence 
model,  investigates  the  effect  on  the  crossing  frequency  of 
different  spectral  functions:  von  Karman,  Kaimal  and 
Teunissen  for  the  x  and  z  directions  with  five  filters: 
ideal  band-pass,  ideal  low-pass,  "quadratic-type",  "sine- 
type"  and  the  Banning.  The  filters  have  a  much  greater 
effect  than  the  spectral  functions.  The  estimated  crossing 
frequency  variation  is  as  much  as  50  percent  among  the 
quadratic-type,  low-pass  and  sine-type  filters.  The  Banning 
filtering  predicts  crossing^  rates  up  to  thirty  times,  and 
the  ideal  one  decade  wide  band-pass  filtering  predicts  rates 
of  eight  times,  higher  than  the  ideal  low-pass  filtering. 
The  variation,  between  the  x  and  z  direction,  is  less  than 
10  percent  for  the  von  Karman  spectrum,  and  over  40  percent 
for  the    Teunissen  one. 
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I.  INTRO DDCTION 

Most  flows  occurring  in  nature  and  in  engineering  appli- 
cations are  turbulent.  In  the  earth's  boundary  layer, 
differential  heating  of  the  atmosphere  produces  pressure 
gradients,  which  are  subsequently  modified  by  the  rotation 
of  the  earth,  causing  a  complex  velocity  field.  In  this 
boundary  layer  the  wind  speed  decreases  as  the  surface  is 
approached  due  to  both  the  frictional  drag  of  the  surface 
and  the  drag  of  all  bodies  protruding  into  the  air  flow. 
These  retarding  forces  are  transmitted  through  the  layer  by 
shear  forces  and  by  the  exchange  of  momentum  due  to  the 
vertical  movement  of  the  air.  The  process  of  momentum 
exchange  between  layers  is  the  mechanism  leads  to  the  gener- 
ation and  decay  of  eddies  which  are  termed  turbulence.  The 
resulting  mixing  of  the  air  produces,  along  all  three 
orthogonal  axes,  fluctuations  in  wind  speed,  commonly  called 
gusts,  which  vary  in  size  in  both  time  and  space. 

The  effects  of  this  atmospheric  turbulence  has  been  of 
continuing  concern  to  the  aircraft  or  structure  designer. 
Typical  turbulence  related  problems  are:  the  effects  of 
turbulence  on  the  fatigue  life  of  the  structure;  the 
performance  of  control  systems  in  turbulence  and  the  deter- 
mination of  ultimate  structural  strength  ret^uired  to  accept 
unsteady  loads  induced  by  turbulence.  In  an  attempt  to  solve 
these  problems,  a  number  of  statistical  models  of  turbu- 
lence, which  endeavour  to  describe  the  turbulence  in  terms 
of  as  few  parameters  as  possible,  have  been  proposed  since 
early  in  this  century.  To  analyze  these  problems,  one  needs 
to  know  the  crossing  frequency,  which  is  the  frequency  with 
which  a  random  function (  here  the  filtered  wind  speed) 
crosses  a  prescribed  value  defined  by  a  certain  averaging 
time,  say  the  hourly  mean. 


The  crossing  frequency  depends  largely  on  whether  or  not 
the  random  function  is  filtered,  the  type  of  filter  and  to  a 
lesser  extent, the  spectral  density  function.  The  effects  of 
type  of  filter  used  and  spectral  density  expression  will  be 
discussed  in  detail  later.  The  spectral  density  functions 
have  been  proposed  by  Dryden  [Ref.  1],  von  Karman  [Eef.  2], 
Kaimal  [Eef.  3]  and  Teunissen  [Eef-  4].  The  Dryden  expres- 
sion makes  mathematical  analysis  simple  but  is  much  less 
accurate  than  the  others.  The  von  Karman  spectral  function 
was  recommended  by  the  Engineering  Sciences  Data  Unit 
(E.S.D.U.)  [Ref.  5]  and  this  expressions  was  originally 
postulated  for  isotropic  turbulence.  The  effect  of  departure 
from  isotropic  turbulence  near  the  ground  is  allowed  for  by 
the  variation  with  height  and  surface  roughness  of  the 
appropriate  variance  cT  and  length  scale  parameter  L|  since 
they  typify  the  intensity  and  size  of  eddies  constituting 
turbulence.  The  Kaimal  expression,  which  is  frequently 
used,  was  obtained  recently  for  the  surface  layer  over  flat, 
uniform,  relatively  featureless  terrain  in  Kansas.  The 
Teunissen  model,  which  is  a  modified  Kaimal  one,  was 
obtained  for  the  generally  rougher  gross  features  of  the 
upstream  terrain. 

In  the  following  discussion,  we  will  review  turbulence 
models  and  spectral  density  functions.  The  spectral  density 
S(n),  in  non-dimensional  form,  is  a  function  of  the  dimen- 
sionless  parameter  nL^/U,  where  n  is  frequency  and  L^*  is  the 
length  scale  of  the  turbulence  and  0  is  the  mean  hourly  wind 
speed.  E.S.D.U.  [ Bef •  5]  gives  empirical  expressions  for 
the  length  scale  obtained  by  analysing  its  collection  of 
world-wide  turbulence  data.  Finally,  for  given  values  of 
nL^/U,  we  will  estimate  the  crossing  frequencies  with  von 
Karman,  Kaimal  and  Teunissen  spectral  expressions. 
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II.  TDRBOLENCE  MODELS 

A.   DISCEETE  COSTS  HODEl 

Powell  and  Connell  [Ref.  6]  defined  gusts  as  constitu- 
ting any  series  of  discrete  velocity-time  events  that  can  be 
defined  from  a  turbulence  time  series  according  to  some 
extrinsic  criterion.  There  are  two  fundamentally  different 
treatments  of  gust  time;  one  is  that  the  gust  time  is  arbi- 
trarily fixed  and  the  other  is  that  the  gust  time  varies. 
Examples  of  this  latter  type  of  gust  definition,  shown  in 
Figure  2.1,  represent  a  sample  of  wind  fluctuation  in  the 
form  of  a  wind-component  time  series. 

The  GUSTO  and  GUST1  models  define  gust  events  in  terms 
of  gust  amplitudes  and  characteristic  times.  Definitions  of 
amplitudes  and  times  differ  for  the  two  models.  The  GUSTO 
model  can  be  completely  specified  in  terms  of  a  positive 
peak  amplitude  +A<,  ,  and  the  time  interval  ♦T,,  between  zero 
crossings  on  either  side.  It  can  also  be  specified  in  terms 
of  a  negative  peak  amplitude  -Aq  and  the  corresponding  time 
interval  -Tq  .  Both  definitions  are  expected  to  yield 
similiar  statistics.  Powell  and  Connell  combine  both  posi- 
tive and  negative  values  into  one  set  for  their  analysis. 

The  GUST1  model  can  be  specified  in  terms  of  a  positive 
amplitude  ♦Aj  which  is  the  peak-to-peak  amplitude  between 
adjacent  minima  and  maxima,  and  a  time  interval  +Ti  between 
the  minima  and  the  maxima;  the  positive  sign  indicates  a 
positive  rate  of  change  in  the  variable.  A  comparable  defi- 
nition can  be  made  in  terms  of  amplitude  -Aj  and  time  -Tj  , 
where  the  negative  sign  indicates  a  negative  rate  of  change 
in  the  variable. 
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GUST  0  Definition 


GUST  1  Definition 


•*T, 


Figure  2.1    Definitions  of  GOSTO  and  G0ST1. 

The  discrete  model  treats  turbulence   as   a  series  of 
isolated  gusts  and  the  basic  assumptions  are 

1)  Atmospheric  turbulence  can  be  modeled  as  a  collection 
of  isolated  gusts  randomly  distributed  along  the 
structure. 

2)  Gusts  have  random  magnitude  but  fixed  shape. 

3)  The  structure  is  a  deterministic   linear  system  with 
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sufficient  damping,  that  each  encounter  with  one  of 
these  gusts  results  in  a  single  significant  response 
peaks. 
This  model  is  typically  used  to  estimate  ultimate 
strength  requirements,  however,  it  is  not  a  satisfactory  one 
for  all  aspects  of  the  design  problem.  For  example,  its 
application  in  areas  such  as  structural  fatigue,  control 
system  performance,  or  even  extreme  responses  which  involve 
lightly  damped  modes  is  highly  questionable.  In  most 
instances  turbulence  does  not  occur  as  discrete  gusts  but  as 
a  continuous  random  disturbance.  Furthermore,  this  model 
neglects  most  of  the  dynamic  characteristics  of  the  vehicle. 
Thus,  although  the  necessity  of  evaluating  the  responses  of 
a  proposed  structure  to  discrete  gusts  is  still  recognized 
eis  an  important  part  of  the  design  procedure,  the  use  of  the 
discrete  model  to  calculate  most  response  statistics  has 
largely  given  way  in  recent  years  to  the  use  of  turbulence 
models  which  attempt  to  take  both  the  dynamic  characteris- 
tics of  the  structure  and  the  continuous  nature  of  turbu- 
lence into  consideration.  This  is  known  as  the  continuous 
or  power  spectral  model. 


B.   CONTINOOOS  MODEL 

The  basic  idea  of  continuous  model  is  that  atmospheric 
turbulence  can  be  represented  by  a  continuous  stochastic 
process  which  acts  as  a  disturbing  influence  on  the  struc- 
ture. The  principal  assumptions  for  this  model  are 

1)  Each  encounter  of  an  structure  with  continuous  atmos- 
pheric turbulence  can  be  modeled  as  a  deterministic 
linear  system  perturbed  by  three  independent  stati- 
onary stochastic  processes,  which  represent  the  long- 
itudinal, lateral  and  vertical  gusts  components. 

2)  The  spectral   density   of  each   random  process  has  a 

13 


form  characterized  by  three  parameters,  the  gust  vel- 
ocity standard  deviation  ((5*)  ,  the  scale  length  (L)  ,  and 
the  averaged  speed (U) . 

3)  Each  of  the  three  gusts  components  is  a  Saussian 
process. 

The  assumptions  of  the  Gaussian  model  make  it  possible 
to  calculate  the  statistics  of  any  structure  response  for 
each  region  of  turbulence  as  functions  of  the  parameters  L 
and  C  of  that  region.  Thus,  the  assumptions  of  structure 
linearity  and  the  Gaussian  nature  of  turbulence  permit  their 
evaluation  with  a  minimum  of  difficulty.  These  results  will 
be  dependent  upon  the  assumed  values  of  L  and  o*  along  with 
the  characteristics  of  the  structure  dynamics. 

The  continuous  model  is  used  primarily  to  evaluate 
response  statistics  for  selected  conditions  in  continuous 
turbulence.  For  this  application  both  the  scale  length  and 
standard  deviation  of  the  turbulence  as  well  as  the  struc- 
ture characteristics  are  fixed  at  values  representative  of 
certain  conditions  which  could  reasonably  be  expected  to 
occur  in  service. 

The  three  assumptions  made  above  need  some  reconsidera- 
tions. The  first  assumption  stated  that  the  turbulence  is 
stationary  and  homogeneous  and  the  dimensions  of  structures 
are  much  smaller  than  the  scale  lengths  of  the  gusts  compo- 
nents. These  conditions  are  not  always  satisfied.  The  second 
assumption  of  the  power  spectral  model  is  valid  at  low  alti- 
tudes, but  is  not  so  certain  at  high  altitudes,  the  problem 
made  being  more  difficult  by  a  lack  of  data  at  the  latter 
altitudes.  The  third  assumption  stated  that  the  turbulence 
is  a  Gaussian  process.  However,  some  experimental  data  indi- 
cate that  the  Gaussian  turbulence  model  significantly  under- 
estimates the  frequency  of  occurence  of  high  gust  velocities 
[Ref.  7].  Furthermore,  peak  gust  velocity  cumulative  prob- 
abilities of   exceedance  observed  in  atmospheric  turbulence 


1U 


iehave  like  exp  (-x)  rather  than  exp{-x^),  as  predicted  by 
the  Gaussian  distribution  [Ref.  8].  Recently  Reeves,  et  al. 
[fief.  9]  have  discussed  in  great  detail  a  non-Gaussian 
model.  On  the  other  hand,  Ref erences  (6)  and  (10)  indicate 
that  filtered  wind  speeds (particularly  with  a  band-pass  5-50 
Hz    filter)         is   reasonably  Gaussian      in  character.  In   the 

present   work,    the   Gaussian   model  is    used   exclusively. 

The  statistical  quantities  of  most  frequent  interest  are 
the  mean,  variance,  probability  density,  autocorrelation 
function,  power  spectral  density  and  crossing  frequency.  The 
power  spectral  shapes  usually  assumed  are  those  proposed  by 
von  Karman,  Kaimal  and  Teunissen.  The  crossing  frequency 
with  these  three  spectral  expressions  will  be  studied  in 
more    detail   later. 
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III.    SPECTRAL    DENSITY    AND    PROBABILITY    DENSITY 

A.       SPECTRAL   DENSITY 

In  the  field  of  mechanical  and  electrical  oscillations, 
a  large  number  of  phenomena  are  governed  by  equations  of  the 
form 


"  3^^    *    =  Tt"    *  "^  =  ^<*' 


where,  m,  c  and  k  are  constants;  f (t)  is  a  given  forcing 
function  and  y (t)  represents  the  response  of  the  system  as  a 
function  of  the  time  t. 

If  the  above  time  dependent  force  f  (t)  is  periodic,  it 
can  be  expressed  in  terms  of  a  Fourier  series;  this  series 
represents  the  sum  of  a  large  number  of  forces  of  different 
amplitudes  and  frequencies.  All  elastic  bodies  have  natural 
frequencies,  and  if  this  force  is  applied  to  an  elastic 
system,  and  the  amplitudes  of  the  components  near  a  natural 
frequency  are  not  very  small,  then  this  force  will  drive  the 
body  into  resonance.  The  body  will  then  go  through  a  very 
large  number  of  stress  cycles  in  a  short  period  of  time, 
possibly  leading  to  fatigue  failure. 

To  have  stationary  properties,  a  random  signal  must  be 
assumed  to  continue  over  an  infinite  time,  and  in  such  a 
case  neither  the  real  nor  the  imaginary  part  of  the  Fourier 
transform  converges  to  a  steady  value.  In  this  case,  we  use 
the  spectral  density  function,  which  has  no  convergence 
difficulty  and  which  is  applicable  to  a  whole  class  of 
similarly  generated  functions. 
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Spectral  functions  of  atmospheric  turbulence  provide 
information  on  the  frequency  distribution  of  the  kinetic 
energy  of  the  various  fluctuating  velocity  components.  Used 
in  conjunction  with  certain  transfer  functions,  they  provide 
information  about  the  dynamic  loading  on,  and  response  of, 
building  and  aircraft  structures  in  the  atmospheric  wind. 

A  considerable  number  of  measurements  of  the  power 
spectra  of  gust  velocities  in  the  near-neutral  atmosphere, 
or  strong  wind  conditions,  at  varying  heights  and  for 
different  terrains  are  available.  The  von  Karman  spectral 
expressions,  which  were  exclusively  used  by  E-S.D.O. 
[Bef.  5],  for  each  velocity  component (x, y  and  z  direction) 
are  given  by 

nSg  (n)  _    4nu  (3.1) 


nSy  (n)  _  4nv(1-»'755.2nv)  (3.2) 


cri  (1  +  283.2nJ.)"'^'^ 


s~2 


nSu;(n)  _  4ntu  (1^-755. 2nw)  (3.3) 

<rj    ~    ( 1  +  283. 2n5)"/^ 

where  ny=  Lun/U,  n^=  Ly.n/U,  n^=  L^n/O  and  L^*  is  integral 
length  scale  which  is  a  function  only  of  the  terrain  rough- 
ness and  height  above  ground  for  the  x  and  y  directions;  for 
the  z  direction  it  is  a  function  of  elevation  only,  c^i  is 
the  variance  of  the  wind  velocity  fluctuations  about  the 
hourly  mean.  These  equations  are  in  common  use  in  the  engi- 
neering literature  and  have  the  advantage  that  the  integral 
length  scale  L^'  are  treated  as  'free'  scaling  parameters 
which  are  chosen  to  match  the  estimated  scales  for  a  partic- 
ular height  and  terrain  type,  while  maintaining  constant 
spectral  shape. 
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Another  commonly  used  spectral  expression  is  the  Kaimal 
one;    £Ref.    3] 

T\Si  (n)       _        0.164  (f/£o)  (3.4) 

<r^*  1+0.  16a(f/f^  /-^^ 

where  nS  (n)  =  logarithmic  power  spectral  density 
(fl        =   variance  of  i 
f    =  reduced  frequency 

fo    =  reduced   frequency   at   the   intercept   of  the 
extrapolated  inertial  subrange  slope  with  the  nS^'  (n)  / 
(5/  =1  line;  or  in  our  notation,  fol   =0.041z/L^'. 
The  Kaimal  expression   is  mere  common  in   the  meteorolo- 
gical literature  and  was  obtained  as   a  best-fit  to  surface- 
layer    measurements   over   uniform,     flat,     relatively 
featureless  terrain. 

Teunissen  [  Ref .  4]  observed  that  von  Karman  model- 
appears  to  be  somewhat  better  than  the  Kaimal  model  at  low 
frequencies.  Tlie  Kaimal  model  largely  underestimates  the 
spectral  content  at  these  frequencies  for  all  compo- 
nents (x,  y, and  2).  On  the  other  hand,  the  Kaimal  model  shape 
appears  to  be  better  than  that  cf  the  von  Karman  model. 

Teunissen  has  proposed  modified  Kaimal  spectral  expres- 
sions, given  by 


(3.5) 
(3.6) 
(3.7) 


nSu  (n) 

= 

10  5f 

».^ 

(0.4U*33f)^'^^ 

nSv  (n) 

= 

17f 

"/ 

(0.38+9.5f/'^3 

nSH^(n) 

= 

2f 

u.^ 

0.a4+5.3f^/^ 
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where  f=nL^*/U  and  u^  is  friction  velocity.  These  expres- 
sions were  derived  from  the  data  which  were  obtained  from  a 
horizontal  array  of  tower-mounted  propeller  anemome- 
ters (z=11in)  during  a  five-hour  period  for  which  the  mean 
wind  direction  was  virtually  perpendicular  to  the  main  span 
of   the   array,    in   the    city   of    Toronto,    Canada. 

As  for  the  reason  for  the  higher  spectral  energies  found 
by  Teunissen  in  the  low-f reguency  region  of  the  spectra  in 
comparison  with  the  Kaimal  values  which  was  measured  in 
Kansas,  this  is  not  entirely  obvious,  apart  from  the  rougher 
terrain   for   Teunissen  experiment.  Teunissen  has    suggested 

that  the  power  spectra  over  a  wide  range  of  terrain  types 
may  possibly  be  represented  by  using  above  expressions  with 
appropriate   values    of  such    'terrain    scaling'    parameters. 

B.       PROBABILITY    DENSITY 

In  general,  a  knowledge  of  the  spectral  density  of  a 
random  process  does  not  enable  us  to  determine  its  prob- 
ability density  or  distribution.  However,  a  common  assump- 
tion which  is  reasonable  in  many  applications  is  that 
atmospheric  turbulence  is  a  'normal'  or  Gaussian  process,  as 
discussed  above,  with  a  probability  density  function  for 
which 


P(i)=    .  i-       exp[-i^/2(n^]  (3.8) 


^iilTL 


where  i=u(t),  v  (t)  or  w(t)  are  wind-speed  fluctuations, 
which  have  a  zero  mean.  Note  that  Equation  (3. 8)  depends  only 
on  the  standard  deviation  of  the  gust  velocity.  Thus,  if  the 
probability  distribution   of  cr     is  known,    the  distribution 
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function  of   the      gust  velocity    can   be      found   by   integration, 
i.e.  , 


P  =      p(i)di  (3.9) 


In  practice,  atmospheric  turbulence  contains  patches  of 
a  significantly  non-Gaussian  nature  (particularly  in  the 
lower  30m)  when  larger  gusts  and  longer  lulls  may  occur  more 
frequently  then  indicated  by  the  Gaussian  distribution. 
Kaimal  et  al.  £Ref.  10]  observed  that  the  assumption  of  a 
stationary  Gaussian  process  has  validity  only  for  the 
filtered  data,    not   for  the   unfiltered   time    series. 
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IV.    CBOSSING    IlEQOENCIES 

The  instantaneous  horizontal  wind  speed  is  considered  to 
he  composed  of  a  mean  longitudinal  component  U  plus  fluctu- 
ating component  v^ ,  where  the  subscript  i  designates  the  x, 
y,  or  z  component  in  a  coordinate  system  with  x  oriented 
along  the  mean  wind  vector  and  z  upwards.  Then  the  wind 
speed   as  a   function   of   time    is    given  by 

7    (t)    =    U   +    v^.  (t)  (4.1) 

The  average  frequency  of  positive-slope  level  crossing 
"Soi ,  with  which  the  wind  speed  fluctuations  exceed  their 
zero  value,  for  a  Gaussian  distribution  of  turbulence,  is 
given  by 

r   r""  /  r"^^  ]^^ 

No^    =       I  n^F(n)S^-  {n)dn  /     /F(n)  S^-  (n)  dn  (4.2) 

^   Aci  ^ci  -' 

where  n^  is  cut-off  frequency,   i  is  u, v,or  w  and  F (n)   is  a 
filter  function. 

The  value  of  Hqi  depends  on  the  mean  wind  speed  U,  L^  and 
turbulence  intensity  a^/U,  which  is  a  measure  of  the  magni- 
tude of  turbulence  fluctuations  and  defined  as  the  ratio  of 
the  standard  deviation  of  the  instantaneous  fluctuating 
velocity  component  to  the  mean  wind  speed  averaged  over 
arbitrary  time,  usually  hourly,  and  depends  only  on  the 
height  z  and  the  terrain  roughness  z©  #  in  a  stable 
atmosphere. 
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A.   FILTER  FUNCTIONS 

It  is  clear  that  the  quality  of  the  crossing  frequency 
estimates  obtained  by  means  of  Equation (4. 2)  will  depend 
strongly  on  the  form  of  the  filter  function  F(n)  and,  hence, 
on  the  selection  of  the  time  function  h{t).  Thus  we  will 
consider  some  examples  of  the  function  h (t)  and  filter  func- 
tion F  (n)  corresponding  to  then.  Let  the  time  function  be 
given  by 


hj^(t)  =  2   sin  (27tt/T)  (4.3) 

T     2Tlt/T 

If   we   take    Fourier   transform,    the   Equation  (4.3)     becomes 

Fj^  (n)    =f    1  when       jnl    <    1/T  (4.4) 

1/2  when       in|    =1/T 

0  when       inj    >1/T. 

Graphs  of  the  functions  h^  (t)'  and  F^  (n)  are  shown  in  Figure 
4.1  and  Figure  4.2  respectively  and  they  represent  an  ideal 
band-pass    filter. 

The    simple    impulse  response  function   h^  (t)       is   the    ordi- 
nary   moving   average: 


1^2^^)    =f  VT        when         -1/2   <t<    +T/2  (4.5) 


'■{ 


0  when  1 t |    >    T/2 


The  corresponding  filter  function  F^  (n)  (F_^  (n)  =H2  (n)"^  )  ,  where 
H^(n)  is  the  Fourier  transform  of  the  function  h2(t),  has 
the   form 


2 
F^  (n)     =    /   simtTn   \  (4.6) 

V      TCTn         J 
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TIME  FUNCTION 


3.0 


TIME(T) 


Figure   4.1        Tine   Function   h^  (t)  ,    li2  f^)  '    i^S  (t)    and   h^(t)  . 

The  function  hj  (t)  and  Fj^  (n)  are  shown  in  Figure  4.1  and 
Figure  4.2  respectively.  The  figure  shows  that  F2  (n)  is  a 
low-f reguency  filter  with  cut-off  freguency  n^  =1/T.  The 
shape  of  this  filter  shows  that  freguencies  above  n^  are 
transmitted  nonunif ormly  and  those  below  n^  are  completely 
suppressed. 

The   impulse    response  function  h,  (t) ,    called  the   smoothed 
moving   average   is: 


1^3  (t)    = 


-M-»-cos(2Tt:/T)  t      when      -T/2   <t<   T/2 
<  T 

0  when        It)    >   T/2 


(4.7) 


23 


FILTER  FUNCTION 
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©■ 

d 
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w 
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LEGEND 
SINE  FILTER 
QUAD."  FILTER 


I.D.B.P.  FILTER 
MANNING  FILTER 


.r.;.. 


/    ••.' 


..\ 


-3.0    -2.5       -2.0       -1.5       -1.0       -0.5         0.0         0.5  1.0  1.5  2.0         2.5 


3.0 


FREQUENCY(1/T) 


Figure   4.2        Filter   Function   F^  (n)  ,    F^  (n)  ,    F3  (n>    and  F.  (n)  . 


The   corresponding   transfer  function   F,  (n)     has  the   form 
Fj  (n)    =  /  simcTn 


(simcTn>  u^ 

TCTn     /    [TL^-{TrTn)^]^ 


(4.8) 


The  advantage   of  F3  (n)    over  F^^  (n)  ,    is  that  it   has  much 
smaller  side  maxima. 

Another  commonly   used  filter  function  is  the  quadratic 
filter  function  given  by 


F^(n)  = 


n^  +  n=^ 


(4.9) 
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One      can   find      the      impulse      function    of      Equation (U. 9)         by 
inverse   Fourier    transform,    i.e., 

=  -^  Ha}  *n'  )^ccs2itntan  + 
■j^  jl^i    +n2)'^*(jsin2]int)dn  {**-  10) 


The   second   term   of   above   equation   is    zero,      since,      it    is   an 
odd    function.      Thus    Equation (4 .  10)    becomes 


h^(t)    =-^/{nc    +n*)  ^  co£27intdn 
Jo 


or 


h^(t)     =-^  Ko  (n^t)  (4.11) 


where  K^ is  the    modified   Bessel    function. 

The  functions  h^  (t)  and  F^  (n)  are  shown  in  Figure  4.1 
and   Figure   4.2. 

The  cut-off  frequencies  are  defined  formally  as  the 
frequencies  at  which  the  filters  reach  the  first  zero.  But, 
when  a  real  filter  is  used,  we  must  find  some  equivalent 
cut-off  frequency,  which  is  defined  as  the  cut-off  frequency 
of  an  equivalent  ideal  filter  whose  area  is  equal  to  that  of 
the  real  filter.  [Bef.  11]  Greenway  [Ref.  12]  and  Powell 
and  Connell  [fief.  6]  have  used  the  value  at  the  half- 
power  {3db)  point  of  the  filter  as  the  cut-off  frequency. 
Vinnichenko  et  al.  [fief.  11]  defined  the  equivalent  cut-off 
frequency   as 

n^    =       r|F(n)|dn    .  (4.12) 
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EFFECTS  OF  CUT-OFF  FREQ. 


10"  16 

CUT-OFF  FREQUENCYjr,^ 


Figure   4.3        Effects  of  Cut-off  Preguencies. 

According  to  the  Equation (4. 12) ,  the  cut-off  frequencies 
depend  on  the  filtering  function-  The  effects  of  cut-off 
frequency  on  the  crossing  frequencies  are  shown  in  Figure 
U. 3,  which  shows  that  the  higher  the  cut-off  frequency  leads 
to   lower  zero  crossing  rates. 
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B-   CAICDLATION  OF  CROSSING  FREQUENCIES 

In  this  section  the  frequency  of  zero  crossings  are 
estimated  using  von  Karman,  Kaimal  and  Teunissen  spectral 
density  functions  with  various  filter  functions  which  we 
have  discussed  above. 

''  •   gsinq  the  von  Karman  Spectral  De_nsity  Function 

Introducing  the  von  Karman  spectral  density  function 
for  the  x-direction  S^ (n)  into  Equation (4. 2)  yields 


'o;j 


Pn^  F  (n)dn        /  P     F(n)  dn 
J(1  +  70.8nif/V  J  (1+70. 8ni 


■,'/2 


f" 


For    numerical  calculation,    set    y=n/nc,  and  D=nc2  /10nc| 
Equation  (4, 13)     transformed   to 


(4.13) 


and 


Nox="ci 


P         y    dy  /  P         dy 

/(UTO.aH^^,  y^)*/y   j(U70.8n/, 


.Vo 


l^f 


(4.14) 


where  Nqx     is  the   dimensionless   frequency      of   positive-slope 
level  crossings    for    the   x   direction    No;(Lu/U.      Equation  (4.  14) 
is      for   the     ideal   band-pass      filter   function      of    D      decades 
filter   width. 

The  upper  integral  of  Eg.  (4.13)  will  not  converge  as 
^C2  goes  to  infinity  because  the  empirical  expression  for 
the  turbulence  spectrum  employed  does  not  account  for  the 
high-frequency  range  where  viscous  dissipation  damps  out 
turbulence  fluctuations.  One  can  achieve  convergence  by  use 
of    a      suitable   filter  function.  Introducing   Equation  (4. 9) 

into    Equation (4. 13)     yields 


ox 


nc      P y'dy  /P  dy 

^(1  +  70.8n/y2)*^'^(Uy*)  /  J(1  +  70.  8ng^  y' P(U 


.^/2 


yM 


(4.15) 
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If  we  take  F^ (n) ,  here  called  the  sine  filter  function, 
instead  of  quadratic  filter,  with  the  cut-off  frequency 
ii.=0.44/T  which  Greenway  [Ref.  12]  has  used,  and  substitute 
into  Equation (4.  6) ,  we  find 

2. 


F^  (n)=sin'^(0.4aT[n/n^)/(0.a4TL  n/n^. ) 


(4.16) 


Introducing  Equation  (4.  16)       into  Ec^uation  (4.  13)       and   setting 
y=0.  44un/nc.    yields 


No-x  =  J^ 


0.  44T1: 


oo 


sin     (y)dy 

[U70.8(nc    /0. 44711)^  y^]   ^ 

[sin^  (Y)/y^]  dy 
[1+70.  8  (nc    /0.44Tt)^  y*]    ^ 


Vx 


(4.17) 


If  we  use  F,  (n) ,  which  is  the  Banning  frequency 
window,  and  equate  cut-off  frequency  n^=0-375/T  as  derived 
by  Vinnichenko  et  al.  [Ref.  11],  then  substitution  into 
Equation  (4.  8)    gives 


F3(n)  = 


siD(0.3757Cn)] 


■n:' 


(4.18) 


2l2- 


0.375itn       /[TL''-(0.375TLn)*] 


Substituting   Equation  (4. 18)     into   Equation  (4. 13)       and   setting 
y=0.375Tin/nc   and     "^^  =n'c/0.  37  5TL  yields 


(y)dy 


yzj% 


(X^-y^)^ 


rCsin^(^ 
o/(1+70.8'7 


(y)  /y^3  <3y 


nci^'f^^  (Tt'-y")^ 


Vz 


(4. 19) 
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TABLE   I 
Crossing  Frequencies (von  Karman   x- direct ion) 


0.  100 

0.  125 

0.  150 

0.  175 

0.220 

0.500 

1.000 

2.000 

5.000 

10.000 

15.000 

22.00  0 

50.000 

75.000 

100.000 

150.000 

175.000 

200.000 

225.000 

250.000 


IDBP 


2UAD 


0.39 

o.as 

0.57 
0.65 
0.81 
1.81 
3.62 
7.23 
18.0 
36.1 
5U.2 
79.6 
180. 
271. 
36  1. 
54  2. 
633. 
723. 
814. 
904. 


72 

38 

09 

85 

72 

93 

35 

94 

931 

84  7 

766 

055 

9212 

38  16 

8425 

7637 

2241 

6853 

1465 

6067 


0. 
0, 
0, 
0, 
0. 


1215 
1116 

1250 
1381 
1604 
0.2737 
0.42  9  8 
0.6765 
1.2365 
1.9555 
2.5585 
3.2989 
5.6928 
7.4555 
9.0288 
11.8268 
13.1054 
14.3243 
11.4  933 
16.6196 


SINE 


0.0865 

0.0955 

0.  1117 

0.1232 

0. 1424 

0.2391 

0.3699 

0.5744 

1.0383 

1.6267 

2.1230 

2.7307 

4.6950 

6.1414 

7.4326 

9.7292 

10.7787 

11.7794 

12.9392 

13.6641 


nc=cut-off   frequency 
IDBP=ideal  bana-cass   filter 
QDAD= "quadratic-type"    filter 
SINE="sine-type"   filter 
HANNING=HANN1nG    filter 
IDLP=ideal   low-pass   filter 


BANNING 


17. 

21 

26 

30 

38. 

87 

174. 

348. 

872. 

1744 

2616 

3836 

3720 


480 

074 

670 

267 

741 

049 

4062 

8093 

0234 

.0481 

.0669 

.9021 

.2461 


The  results  of  numerical  calculations  of  Eqs.(4.14),  (4.15), 
(4.17)  and  (4.19)  are  given  in  table  I.  The  plots  for  above 
equations  are  shown  in  Figure  4.4.  The  plots  showing  how 
well  these  crossing  frequencies  co-incide  depend  on  the 
length  scale,  mean  wind  speed  and  filter  functions. 

The  von  Karman  spectral  density  function  for  the 
z-direction  is  given  by  Equation (3. 3) .  The  equation  of  the 
crossing  frequency  for  the  ideal  band-pass  filter  function 
can  be  obtained  by  substituting  Equation (3. 3)  into 
Equation  (4.2)  and  is  given  by 
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No2=ac» 


I 


(U755.2DU  y^)  y^dy 


'^   2 


»     (1-K283.  2iic>  y   ) 


"/6 


M-t-755.2Dc'y')dy 
(1+283.  2nc?y')"/6 


Vz 


(4.20) 


VON  KARMAN(X-DIRECTION) 


10"        itf 
CUT-OFF  FREQUENCY(nc) 


Figure  4-4    Crossing  Frequencies  with  von  Karman 
Spectral  Expression  for  the  i-direction. 
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Using  the  same  procedure  as  for  the  x-direction,  the 
z-direction  crossing  freguencj  for  the  quadratic  filter 
function  is   found   to   be 


Noz=nc 


00 

r    (1-»-755.2ncV)  y^dy 
-1(1+283.  2 nc'y^f^  (1+y^ 


.Vz 


U 


rT  2  it2 


(H755.2nc''y^)  dy 
( 1+283. 2ncy^)^(1+y^) 


(4.21) 


and,         for      the      sine        filter      function,      if      one        defines 
n^/0.44TC=^g  the   crossing      frequency    for      the    z-direction 

becomes 


r      rCl-^755.2(ff,  yf]sin" 
No2=^J    ^o        (1  +  283. 2n/y^)"-^^ 


(y)'iy 


■  i  (1  + 


2(ncynCsin^  (y)/l^Idy 


(1  +  283. 2»f^2y2)"/6 


iVl 


(4.  22) 


Introducing  filter   function   F^     (n)       into   Equation  (4. 20)       and 
setting   y=0.375Ttn/n^    and   n^. /O.  3  75Tt='tc    yields 


f    rCl^755.2(?cy)^]sin^  (y)  dy 
Noz  =7c  i  "i>      (1+28  3.  2^/y^)"^^(Tl^-y^)^ 

r[1+755.2(n/y'^3[sin^(y)/y^]dy 
,  J       (1+283.  2^,^y2)"''V*-y=')^ 


V. 


(4.23) 
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TABLE   II 

Crossing   Freguenc 

ies(voD  Karman   z-direction)                   1 

nc 

I^BP 

2UAD 

SINE 

BANNING 

0.  100 

0.3812 

0.  121 

0.0816 

17.4448 

0.  125 

0.4686 

0.  136 

0.0934 

21.8044 

0.  150 

0.5567 

0.  150 

0.1044 

26.  1641 

0.  175 

0.6454 

0.163 

0.1144 

30.5238 

0.220 

0.8059 

0.  185 

0. 1320 

38.3716 

0.500 

1.8137 

0.297 

0.2188 

87.2036 

1.000 

3.6207 

0.452 

0.3375 

174.4053 

2.000 

7.2380 

0.696 

0.5238 

348.8093 

5-000 

18.0925 

1.248 

0.9458 

872.0229 

10.000 

36.1844 

1.951 

1.4870 

1744.0461 

15.000 

54.2764 

2.535 

1.9407 

2616.0693 

22.00  0 

79.6054 

3.246 

2.4978 

3836.8989 

50.000 

180.9214 

5.493 

4.2994 

8720.2266 

75.000 

271.  3823 

7.  127 

5.6256 



100.000 

361.84  30 

8.627 

6.8094 



150.000 

542.7634 

11.129 

8.9150 



175.000 

633.2231 

12.330 

9.8773 



200.000 

723.6838 

13.474 

10.7948 



225.000 

814. 1440 

14.572 

11.6748 



250.000 

904.6057 

15.63C 

12.5228 

._    _ - 

The  results  of  numerical  evaluations  of  Eqs.(4.20), 
(4.21)  ,(4.22)  and  (4.23)  are  listed  in  the  Table  II.  The 
plots  for  above  equations  are  shown  in  Figure  4.5. 

For  comparison,  the  crossing  frequencies  with  von 
Karman  spectral  expression  for  the  x  and  z  directions  are 
plotted  in  Figure  4.6.  The  crossing  frequency  for  the  z 
direction  is  less  than  the  x  direction  except  ideal  band- 
pass filter  case,  in  which  the  frequencies  are  identical  for 
both  directions. 

The  crossing  frequency  plots  presented  in  'Figures 
4.4  -  4.6  offer  a  detailed  view  of  how  the  crossing 
frequency  depends  on  the  filter  function.  Thus,  in  designing 
a  system,  the  type  of  the  filter  has  to  be  chosen  carefully 
to  match  the  behavior  of  the  particular  system. 
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VON  KARMAN(Z-DIRECTION) 


10"  lO" 

CUT-OFF  FREQUENCV{nc) 


Figure  4.5        Crossing   Frequencies   with   von   Karman 
Spectral  Expression   for   the  z-direct ion. 


2-      UsiHSl  th.§  Kaimal  Spectral    Densitx   Function 

From  Equation  (3. 4)  ,    the  asymptotic   spectral   behavior 
can    be   written  as 

■V3 


nS,-  (n)   ^r        (f /f ^ )  '='       ;   f»fe 
cr;-*  I  0.164  (f/fo)     ;    f«fo 


(4.24a) 
(4.24b) 
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VON  KARMAN(X  AND  Z  DIRECTION) 

— Ifl 


CUT-OFF  FREQUENCY(nc) 


Figure   U.  6        Crossing   Freguencies  with  von   Karman 
Spectral   Expression  for    the  x  and  z-direction. 


The   expression  for   the    low   frequency   behavior   of   the 
spectrum   is  given 


^^    =  (  4U72)f 


(4.25) 


where   L^*    represents   the   length   scale.      Equating 
Equation  (4.25)  and  Equation  (4.  24b)  becomes 


L^-  =  0.04l2/(fo)„c 


(4.26) 
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lor   the   x-direction    only 
Lu  =   O.Oaiz/foa 
or  foy=   0.041z/Lu 

Substituting      Equation (4. 27)    into    (3.U)    yields 


(4.27) 


nSu  (n)     _        0.164(f/Q.  041z/Lu) 


or 


nSu(n) 


1+0.164(f/0. 041z/Lu) 
4n 


1+33. 641n 


;r^/3 


(4.  28) 


where      n   =      nLy/D.  Then      the  crossin^g      frequency    for      the 

X-direction   becomes 


No*  = 


r- 


00 


F(n)n*dn         /    f     F  (n)  dn 


(U33.641n'/3)    /    J(l+33. 641n^''3) 


Vz 


(4.29) 


Substituting  the        quadratic  filter 

Equation  (4.29)    yields 


function  into 


N 


OX 


/- 


n^dn 


iVi 


NO 


(U  33.  64  In ^-^3)  (n^^  »n^) 
dn 


00 


(U33.641n%)  (eJ    *n^) 


(4.30) 


Eeplacing  n/n^  by  y  in  Equation  (4.30)  leads  to 


No;c  =^c  { 


-^[l-«-33.  641  (n^  y)^3]  (Uy^) 


r  ay 

,J[l  +  33.641(ney)^^^](1  + 


y*) 


(4.31) 
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Setting   y=e        in   Equation (4.  31)    gives 


Noj(    =^c 


oo 


i^^  dx 


[h 33.  641  (nce^)^/3]  (Ue^^) 
J  [1+33.64  1  (n<.e'^)^/3](Ue^'^) 


y/2 


(4.32) 


If        we        introduce        the  sine        filter        function        into 

Equation  (4. 29)  ,    the    crossing   frequency    becomes 


ox 


sin^  (0  .  44TC  n/nc  )  n^  dn 
(1+33.641n^/3)  (0.44Tt/nc)^  n^ 


sin^  {0.4  4Ttn/nc)  dn 


(1+33.641n^/3)   (0.4  4Tt/ne)^  n^ 


Vl 


(4.33) 


Substituting   0.44itn/nc   for   y      in  Equation  (4. 33)  ,      and   define 
n,./0.441l=>]c#    yields 


No^  =^c 


if 


sin    (y)  dy 


nV^ 


4     U33.641  (yfcY)^^^ 
[sin^  (y)  /y^]  dy 
1  +  33.641  (^j.y)^/3 


(4.34) 


For   the      Banning   filter,      substituting      Equation (4 . 1 8)       into 
Equation  (4.28)  ,      setting  y=0.375Tcn/nc  and  define 

nc/0.375Tl=»fc  ,    yields 
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No^=^c 


to 


sin^  (y)<iy 


w 


[1^33.641  ( ^,  y )^/3]    (Tl^-y^)^ 

dy 


(4.35) 


r""  Csin^(y)/y^3dy 
.   -/[1-^33.641  (Pf^yf/^] 


(Tt'-y^f 


TABLE    III 

Crossing   Freq 

uencies (Kaimal  x-direction) 

He 

IDBP 

QUAD 

SINE 

BANNING 

0.  100 

0.4056 

0.134 

0.0868 

17.4532 

0.  125 

0.4938 

0.  152 

0. 1001 

21.8126 

0.  150 

0.5819 

0.  169 

0.1126 

26. 1720 

0.  175 

0.6702 

0.  184 

0. 1243 

30.5316 

0.220 

0.8295 

0.210 

0.  1441 

38.3790 

0.500 

1.8305 

0.340 

0.2445 

87.2088 

1.000 

3.6321 

0.518 

0.3811 

174.4094 

2.000 

7.2456 

0.799 

0.5943 

348.8120 

5.000 

18.0968 

1.434 

1.0744 

872.0224 

10.000 

36.1872 

2.249 

1.6879 

1744.0486 

15.000 

54.2786 

2.931 

2.1987 

2616.0667 

22.00  0 

79.6070 

3.769 

2.8249 

3836.9033 

50.000 

180.9221 

6.471 

4.8291 

8720,2461 

75.000 

271.3823 

8.464 

6.2923 



100.000 

361.8423 

10.243 

7.5860 



150.000 

542.7634 

13.324 

9.8619 



175.000 

633.2239 

14.733 

10.8962 



200.000 

723.6846 

16.102 

11.9076 



225.000 

814. 1450 

17.414 

12.8333 



250.000 

904.6052 

18.678 

13.7017 

The  results  of  numerical  calculations  of  Eqs. (4.29), 
(4.32),  (4.34)  and  (4.35)  are  given  in  Table  III  and  the 
plots  for  above   equations  are  shown    in   Figure   4.7. 

It  is  clear  from  the  Eg-  (3.4)  and  Eg.  (4.26)  that, 
for  the  Kaimal  spectral  density  function,  the  difference  in 
crossing  frequencies  between  the  x  and  z  directions  is  only 
in  the  length  scale  L^ .  Thus  the  plots  for  the  z-direction 
are    exactly  same  as    the   x-direction. 
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Figure  4.7   Crossing  Frequencies  with  Kaimal 
Spectral  Ezpression  for  the  x-direction. 


3.   Using  the  Tennissen  Spectral  Density  Function 

The  modified  Kaimal  spectral  density  expressions, 
vhich  were  derived  by  Teunissen,  are  given  in  Egs.  (3.5), 
(3.6)  and  (3-7).  These  equations  were  derived  by  adding  a 
scaling  parameter  to  each  of  the  Kaimal  spectral  equations 
and  evaluating  these  parameters  from  the  observed  data. 
These  equations   are  function  of  friction  velocity  u^  .   In 
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order  to  estimate  the  crossing  fregaency  as  we  have  done 
earlier,  we  need  to  change  that  equations  as  a  function  of 
variance.  In  E-S.D.  U.  [Ref.  5],  the  relations  between  6'^ 
and  u^  are  defined  as 

0.ac52/u^=(<3;-/U)log  (z/zj  (4.36) 

Teunissen  derived   the  average   values   of  oj/u^,    given   by 

ot^/Ujt    =      2.  84  (4.37a) 

cf^/n^   =      1.92  (4.37b) 

<5'^/M^    =1.27  (4.37c) 

From    Equation (4. 36)     and   Equation (4. 37) ,  one   can    get 

uj    (u-component)    =    0.124  (T'u  (4.38a) 

u^    (v-component)    =    0.271    cfj'  (4.3  8b) 

u^    (w-component)    =    0.620   ctJ"  (4.38c) 

and  substituting  Eq.  (4.38)  into  Eqs-(3.5),  (3.6)  and  (3.7) 
yields 

nSu(n)  13n  (4.39) 


cj^  (0.44  +  33n)^'^3 


nSv(n)      _  4.61n  (4.40) 

CJV^  (0.38+9.5n)^''3 


nS»^(n)      _  1.24n  (4.41) 

cr^^  0.44+5. 3n^'^3 


Using  methods  similar  to  the  above,  the  crossing  frequency 
for  the  x-direction,  with  the  ideal  band  pass  filter, 
becomes 
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'ox 


vMta      2 


J    (0.4U+33n 
*»  net 


^nc2 


dn 


'    nci 


V2 


or 


N 


ox    ~     ^Cl 


J  (0.44O3nc,  yf7J(0.44O3nc,yf^ 


(4.42) 


If      one     uses      the      sine     filter     function,        with        cut-off 
frequency   n^^=0-44/T,    the   crossing   frequency    becomes 


Nox  = 


Jn(0 


r- 

J  (0 


sin^  (0.  44Ttn/nfc)  dn 
(0.44+33n)^''5    (0.44Ttn/ncf 


V. 


sin^  (0.  44un/Dc  )  n^  dn  ^  ^ 

4  (0.44-»-33n)^''3    (0  .44nn/nc  )^ 


(4.43) 


Setting   y=0.44Ttn/nc    into  Eg.  (4.43),      and   define    n^ /0. 44Tt=i^c  # 
yields 


r?*       -.-2 


^     r    r      sin^(y)dy  /  fCsin^  (y) /y^]  dy 

Nox='^ciJ       (0-'*'++33?j['^y)^/3/j^   (0.4  4*339i^y)^3 


^V2 


(4. 44) 


For   the   quadratic   filter   function,   IT^^   becomes 


^ox   =    ^c 


Jo(0 


y^dy 


(0.44+33irc  y)^    (1+y^) 


l    J(0 


dy 


(0.44*33nc  y}^/^  (1+y^) 


V, 


(4.45) 
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TABLE   IV 

Crossing   Frequencies (Teunissen   x-direction) 

lie 

IDBP 

^UAD 

SINE 

l^kl 

0.100 

0.3722 

0.060 

0.0425 

0.0334 

0.  125 

0.4628 

0.06S 

0.0494 

0.0394 

0.150 

0.5534 

0.078 

0.0559 

0.0449 

0.  175 

0.6U39 

0.086 

0.0620 

0.0502 

0.220 

0.8068 

0.099 

0.0724 

0.0590 

0.500 

1.8201 

0.16S 

0.1251 

0.  1044 

1.000 

3. 6294 

0.266 

0. 1933 

0.  1672 

2.000 

7.2478 

0.420 

0.3142 

0.2664 

5.000 

18. 1031 

0.768 

0.5754 

0.4915 

10.000 

36.1952 

1.213 

0.9083 

0.7801 

15.000 

54.2873 

1.581 

1.  1846 

1.0219 

22.000 

79.6162 

2.039 

1.5228 

1.3188 

50.000 

180.9321 

3.521 

2.6002 

2.2786 

75.000 

271.3931 

4.612 

3.3801 

2.9851 

100.000 

361.8535 

5.586 

4.0689 

3.6157 

150.000 

542.7749 

7.266 

5.2816 

4.7372 

175.000 

6  33.2354 

8.053 

5.8554 

5.2496 

200.000 

723.6  956 

8.804 

6.2870 

5.7380 

225.000 

814.1565 

9.525 

6.8032 

6.2065 

250.000 

904.6169 

10.219 

7.3011 

6.6579 

The  results  of  numerical  calculations  of  Egs.  (4.42)  , 
(4.44)  and  (4.45)  are  given  in  Table  IV  and  the  plots  for 
these  equations    are    shown   in   Figure    4.8. 

The  modified  Kaimal  spectral  expression  for  the 
z-direction  is  given  by  Eg.  (4.4  1).  If  we  put  ideal  band  pass 
filter  function  and  spectral  expression  for  the  z-directon 
into    Eg.  (4.2)    it    yields 


oz 


/'lea                        /  fj^cz 
n   dn  /  I  dn 

0.44+5.  3n^/5  /  J   0.44+5.3n 
T\ci  /    yic, 


V/z 


^/3 


41 


or 


No£  =    ^ci 


nc2  /  J^c2 

r        y^dy  /  r      dy ^^ 

J  0.44+5.3  (nc,  yfV  ^   0.44+5.3  (He,  yf^ 


(4.46) 


o 
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+     TIDLPX 
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10"  10' 

CUT-OFF  FREQUENCY(nc) 


1         I      I     I    r 


Figure  U.8   Crossing  Fregnencies  using  Teunissen 
Spectral  Expression  for  the  x-direction- 
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If  one  introduces  the  sine  filter  function,  the  equation  for 
the  crossing  frequency  becomes 


'oz 


Jo  (0 


sin^  (O.ttUTin/Pc)  n^dn 
(0.44  +  5.  3nV5)  (0.4  4Un/ncf 


.V, 


sin^  (0.  44iin/nc)  dn 

(0.44+5.3n^/3)  (0.4  4Ttn/ncf  ^ 


(4.47) 


or 


'0£ 


=  \ 


1 
I 


CO 


sin    (y)dy 


'o  0.44+5.3  (He  y)^^^ 
[sin^  (y)  /y^]dy 
0.44+5.3  (rf^y)^''^ 


Vi 


(4.48) 


If   one   introduces   the  quadratic     filter   function,       the    equa- 
tion   for  the  crossing  frequency  becomes 


OH 


=    ^c< 


i,[0.44+5.3(nc 


y)^/3]     (Uy2j 


r^ dy 

,  J[0.44+5,3(n^y)^/3j  p  + 


y^) 


Vi 


(4.49) 


The  results  of  numerical  calculations  of  Eqs.  (4.46), 
(4.48)  and  (4.49)  are  given  in  Table  V.   The  plots  for  above 
equations  are  shown  in  Figure  4.9. 
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TABLE   V 

Crossing   Frequencies (Teunissen  z-direction) 

ic 

IDBP 

QMU 

SINE 

IDLP 

0.  100 

0.4454 

0.178 

0.1098 

0.0559 

0.  125 

0.5367 

0.201 

0.1258 

0.0691 

0.150 

0.6  26  3 

0.222 

0.  1408 

0.0818 

0.  175 

0.7151 

0.241 

0.1551 

0.0  94  1 

0.220 

0.8740 

0.274 

0. 1794 

0. 1153 

0.500 

1.8651 

0.43£ 

0.3047 

0.2268 

1.000 

3.6560 

0.660 

0.4758 

0.3783 

2.000 

7.2610 

1.007 

0.7414 

0.6104 

5.000 

18. 1053 

1.789 

1.3353 

1.  1232 

10.000 

36. 1926 

2.791 

2.0917 

1.7725 

15.000 

54.2827 

3.629 

2.7242 

2.3144 

22.000 

79.6102 

4.658 

3.5000 

2.9785 

50.000 

180.9  239 

7.973 

5.9955 

5.1216 

75-000 

271.3835 

10.410 

7.8197 

6.6982 

100.000 

361.8437 

12.591 

9.4420 

8. 1050 

150.000 

542.7644 

16.4  70 

12.3095 

10.6064 

175.000 

633.2249 

18.242 

13.6158 

11.7493 

200.000 

723.6853 

19.9  32 

14.86  06 

12.8388 

225.000 

814. 1458 

2  1 .  5  53 

16.0474 

13.8837 

250.000 

904.6062 

23.  115 

17. 1816 

14.8905 

For  comparison,  the  crossing  frequencies  with  the 
von  Karman,  Kaimal  and  Teunissen  spectral  expressions  for 
the  x-direction  with  various  filter  functions  are  plotted  in 
Figure  4.10.  As  shown  in  Figure  4.10,  the  crossing 
frequencies  predicted  by  the  Kaimal  spectral  expression  are 
higher   than  for    the   others.  In   Figure    4.11,      the   crossing 

frequencies,  with  all  three  spectral  expressions  and  the 
z-direction,  are  compared.  The  figure  characteristics  are 
similar  to  those  for  the  x-direction,  except  that  the 
crossing  frequencies   are   lower. 
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TEUNISSEN(Z-DIRECTION) 


-1 — l—TT 
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CUT-OFF  FREQUENCY(ffc) 


Figure   U.9        Crossing    Freguencies   with   Teunissen 
Spectral   Density   Function  for   the  z-direction- 


C.       EXAMPLE    CALCDLATIOHS 

In  the  previous  section  of  this  study,  we  derived  equa- 
tions for  the  crossing  frequency.  The  purpose  of  this 
section  is  to  illustrate  the  application  of  these  equations 
to  a  given  turbulence  condition,  and  to  estimate  the 
crossing   frequencies  for   a   low -pass    filter. 
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VON  KARMAN,  Kf  IMAL&  TEUNISSEN(X) 


1 — I — I  1 1 1 1 


T — I — r-TT 


CUT-OFF  FREQUENCY(nc) 


Figure  4.10        Cosparison   of    tozi  Karman.    Kalmal   and   Teunissen 
Spectral   Density  Function  for  the  x-direction. 

The  dimensionless  crossing  frequency,  ^ox  *  is  a  function 
of  n^  (=nl^  /U) .  The  length  scale  parameter  L^  and  surface 
roughness  parameter  Zf,  are  given  by  E-S.D.D.  [Ref.  5].  For 
convenience  these  data  are  included  here  as  Figure  4.12  for 
length  scale  parameters,  and  Figure  4.13  for  the  surface 
roughness  parameters. 
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VON  KARMAN,  KAIMAL&c  TEUNISSEN(Z) 


CUT-OFF  FREQUENCY(nc) 


Figure   4.11        Comparison  of   von  Karman,    Kaimal  and  Teunissen 
Spectral   Density   Function  for   the   z-direction. 


Suppose   we    have    a  site   for    which 
z   =  20m 
nc.=  50Hz 

D  =      25m/sec  and   the      turbulence   condition     is  that 
for    an   off-sea      wind    in  a   coastal   area. 

With  this      choice  of   turbulence   condition,         the   surface 
roughness   parameter   z^   can   be   found      in  Figure    4. 13,      and  is 
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Zo=0.001.      The    length     scale    parameters  are   given      in   Figure 
4.12,    and   are    given    as   L^=105,    1^=36    and   Lu;=6.5. 

The  next  step  is  calculation  of  the  dimensionless 
cut-off   frequencies,    i.e., 

ncu=   nLu/D    =  210  (^.50) 

ntv=   nLv/U    =     72 
n^jy=   nL„;/0    =      13 

The  crossing  frequencies  are  calculated  by  the  use  of 
the  given  crossing  frequency  vs  cut-off  frequency  graphs, 
with  dimensionless  cut-off  frequencies  which  are  given  in 
Eq.  (4.50). 

The  final  step  is  to  convert  the  dimensionless  crossing 
frequency  to  a  dimensional  frequency  by  multiplying  the 
factor   U/L^ ,    i.e.. 

No*  =    Nox  D/Iu  C*.  51) 

NoE  =    Noz   D/Iw 

'  The    results    of   calculation    cf   the   crossing   frequency   for 
a  given   turbulence   condition   are  given   in    Table    VI. 
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TABLE    VI 

Example 

Calculations 

of    the  Crossing   Freg 

uency 

Z2Ji   Karman  ( 

x-direction) 

REMARK 

IIEF 

IDIP                 2UAD 

SINE 

Nox 

Nox      • 

75.0 

17.857 

(ncu=210. 

10,2                 14.6 
2.428              3.476 
0^       U/Lu=0.238) 

12.0 
2.857 

von   Karman ( 

z-direction) 

REMARK 

I  EBP 

IDIP                 2UAD 

SINE 

Noz 
Noz 

4.70 
18.076 

1.50                2.40 
5.769              9.231 
U/Lw=3.846) 

1.80 
6.923 

KaimaKx   an 

d   z-direction) 

REMARK 

IIEP 

IDLP                S^AD 

SINE 

Nox 
Nox 
Noz 
Noz 

75.0 
17.857 
4,70 
18.076 

10.9                 16.50 
2.595              3.928 
1.70                2.65 
6.538              10.192 

12.2 
2.905 
2.05 
7.885 

Teunissen{x 

^H^  z-direction) 

REMARK 

lEBP 

IDIP                QUAD 

SINE 

Nox 
Nox 
Koz 
Noz 

75.0 
17.857 
4,70 
18.076 

6.C8                9.30 
1.448              2.214 
2.10               .3.30 
8.076              12.692 

6.50 
1.548 
2.40 
9.230 

49 


4    5     6 


2    m      (effective   height) 
Vcilues  of   Length   Scale   Parameter  L; 


Figure  U.12 

(Reproduced  from  ESdo  74031) 
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Figure    U.13        Values  of   the   Surface   Roughness   Paraneter   z^ 
(Reproduced  from   ESDO   740J1)- 
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V.    CONCIOSIONS 

This  study  has  estimated  the  crossing  frequency  based  on 
the  three  spectral  density  functions,  which  assumes  turbu- 
lence to  be  a  homogeneous,  stationary  Gaussian  process,  with 
five    different    filter  functions  for    the    x   and  z    directions. 

This  estimation  has  shown,  for  the  x  direction,  that  the 
Kaimal  spectral  function  gives  about  one  percent  higher 
crossing  frequency  than  the  von  Karman,  and  forty  percent 
higher  than  the  Teunissen  one,  at  n't  =  100  for  quadratic-type 
filtering;  for  the  ideal  low-pass  and  sine-type  filters,  the 
difference  between  the  Kaimal  and  von  Karman  spectral  func- 
tion is  only  a  few  percent,  hut  for  the  Teunissen  one  is 
more    than  forty    percent   than   for  the   other    spectrums. 

For  the  z  direction,  the  Teunissen  spectral  function 
gives  about  ten  percent  higher  crossing  frequency  than  the 
Kaimal  one  and  about  twenty  percent  higher  than  the  von 
Karman  one,  at  nc=100  for  the  quadratic-type,  sine-type  and 
low-pass  filtering.  However,  this  difference  becomes  to  zero 
as   non-dimensional   cut-off   frequency    goes    to  zero. 

The  crossing  frequencies  for  band-pass  and  Banning 
filtering  are  only  a  fraction  of  a  percent  different  for  the 
three   spectral   functions. 

The  present  study  shows  that  Manning  and  band-pass 
filtering  predicts  higher  crossing  frequencies  than  low-pass 
filtering.  This  difference  is  greatly  amplified  with  the 
higher   frequency. 

If  some  site  statistics  are  available,  as  for  the 
example  calculation,  then  the  number  of  crossings  that  the 
wind-speed  fluctuations  exceed  their  zero  value  per  a 
certain  period  can  be  computed  and  an  estimation  of  fatigue 
life    can   be   obtained. 
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As  a  consequence  of  the  results  studied  herein,  it  is 
suggested  that,  in  future  work,  careful  consideration  should 
he  given  to  the  type  of  filtering  used  to  model  structural 
behavior  and  determine  more  accurately  the  fatigue  life  of  a 
system. 
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